Thermodynamics of trapped interacting bosons in one dimension 
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On the basis of Bethe ansatz solution of bosons with (5-function interaction in a one-dimensional 
potential well, the thermodynamics equilibrium of the system in finite temperature is studied by 
using the strategy of Yang and Yang. The thermodynamics quantities, such as specific heat etc. are 
obtained for the cases of strong coupling limit and weak coupling limit respectively. 



I. INTRODUCTION 



There has been much interest recently in the study of strongly correlated electronic systems in one dimension. 
This is because not only various non-perturbative methods |]i|~|l0| are applicable to models in one dimension but also 
several photo-emission experiments for one-dimensional compounds of alkali-metal cupper oxides ||ll|-[l3|| are carried 
out. Bosons with (5-function interaction in one dimension is a simple but interesting model. It was solved in ref. 



[|l4| under the periodic boundary condition, in ref. [|15 16 1 under the boundary condition of potential well of infinite 
depth and in ref. ||l^ under that of potential well of finite depth. A strategy for studying the thermodynamics of 
exactly solvable models was proposed in ref. Q when discussing the solution of ref. 0. In the present paper, using 
the strategy of ref. ||l8| we study the thermodynamics on the basis of the Bethe ansatz solution of ref. |17|. After 



recall the model Hamiltonian and the Bethe-ansatz equation we study the thermal equilibrium in Section II. The 
formal expressions of free energy and pressure are obtained. In Section III we consider the strong coupling limit and 
obtain the quasi-momentum distribution and specific heat explicitly. In Section IV the case of week coupling limit is 
discussed extensively. 

II. THERMAL EQUILIBRIUM AT FINITE TEMPERATURE 

The Hamiltonian of bosons in a one-dimensional potential well of finite depth with ^-interaction reads 

^ = -Ear^ + E^(^') + 2c E -5(2;. -X,), (1) 

eq. (Q) is the first quantization version of Gross-Pitaevski equation, which was widely used to study the Bose- 
Einstein condensation [2C[| in recent years. The model Hamiltonian (|l|) was solved by means of Bethe-ansatz method 
1^ . The logarithm of the Bethe ansatz equation reads 



where 



— /, =fc, + -sin-i(^) + iy 



tan {— ) + tan {— ) 



(3) 



where the integers Ij play the role of the quantum numbers. Equation (H) is the secular equation to determine the 
spectrum. Moreover the transcendental equations (|^) is difficult to solve directly. Now, we consider the problem 
in thermodynamic limit: iV 3> 1 and L 3> 1 with a fixed concentration D = N/ L by introducing a "smooth" 
positive-defined density p{k) describing the distribution of roots and holes 
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Treating kj as continuous variable k and differentiating eq(P) with respect to fc, we get an integral equation 



2TTp{k) = 1 



dk'p{k' 



L 



c2 + (fc-fc')2 C^ + {k + k')' 

c c 



C2 + (fc - /l„)2 c2 + (/C + /l„)2 



(4) 



where we have used the replacement 



N 



in the thermodynamics limit. The summation in the right hand side runs over "holes" (including kj) which can be 
written formally as an integral J ph{k)f{k)dk with Ph{k) — i^/ L)'^j^S{k — hm)- Furthermore, eq. (Q) is write out 
as an integral equation for the density of holes ph{k) and the density of roots p{k) — p{k) — ph{k), 



2Tr{p + ph) = 1 + 



dk'p{k')[ 



c^ + {k-k'f c^ + {k + k'Y 



(5) 



As the Bethe ansatz solution is obtained for the case of bounded states (i.e. IvciKj ~ 0), so the range of the 
integration is [— Vq, Vq]. In terms of the distribution function of roots, we can write out the energy per particle 



E/N = / p{k)k^dk, 



where 



D = N/L = J p{k)dk. 
On the basis of the strategy of ref . ||l8| , the total entropy of the system is obtained 

S/N = D^^ I [{p + ph)\n{p + ph) - pin p - ph\nph]dk. 



(6) 



(7) 



(8) 



where the Boltzmann constant is put to unit. 

In the thermal equilibrium, the free energy Q = [E — TS — j-iN) should be in minimmn. Writing 



Ph{k) 
p{k) 



= exp[e(fc)/r], 



we obtain from the minimizing requirements 5^ = the following equations. 



e{k) = -p. + e-^ /[ 



27r 



Eq. (P) is readily written as 



27r(l + e^/'^)p(fc) = 1 + y 



c^^{k-k'f c^ + {k + k'f 



dk'p{k')[ 



]\n{l + e-<''"^/'^)dk'. 



L^/V^~k^ 



c'^ + {k-k'f c^ + {k + k'f 



(9) 



(10) 



(11) 



Principally, eq. (10) can be solved by iteration and then eq. (||) can be a Fredholm type equation for p{k). 

We would like to mention some points about the parameter /i. If minimizing the Helmholtz free energy F = E — TS 
under the condition that the concentration Z? in (|7|) is a constant, one will have a Lagrangian multiplier. The multiplier 
function is just the same as the chemical potential ^ when considering a grand assemble. So both procedures are 
equivalent. 

Multiplying eq. dlO) with pD^^ and integrating over fc, we obtain 



P = 



I - e)pdk [27r(p + p„) - 1 ~ ^ 



ln(l + e-"/'^)dfc. 



(12) 
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The argument k of p and e are always omitted in our notions as long as it does not bring about confusions. With the 
help of eq. the entropy (0) is rewritten as 



N 



= D-^ J [{p + ph) ln(l + e-'/^) + pe/T]dk. (13) 



The Helmholtz free energy per particle is 



^ 1 /(fc^ - e)pdk l[ip + Ph) HI + e-^/^)dk. (14) 



Comparison of eq. (12) and eq. (|lj) gives rise to 

f.,«-_y[i+__=^ii„,i+.-./^)*^ (15) 

Thus the free energy will be obtained once the e(fc) is solved from eq. (p^). As in thermodynamics F = —PL + pN, the 
pressure is P = —{dF/dL)T- It was shown that if e, p are implicit functions of some thermodynamic quantities 
X (such as T, L), the derivative of eq. (15) with respect to x is the same as the partial derivative of eq. ( p^ ) with 



respect to the explicit variable x. Then it is easy to calculate the pressure in terms of the e, namely, 

P-£/[l + i^^]ln(l + e-/-)... (16) 

It is formally similiar to the result of ref. except one more term arising from the boundary condition is involved. 
Like- wisely, the entropy is 

The other thermal quantities such as specific heat etc. are also obtainable if one is able to solve e{k) from eq. (|l^). 
We will calculated them in some special cases in the following sections. 

III. THE STRONG COUPLING LIMIT 

It is difficult to obtain an explicit expression of e(fc) from eq. ( p^ ) generally. However, in some special cases, we are 
able to obtain some plausible results. In strong coupling limit c ^ Vb, eq(pO|) becomes 

e(fc) = -p + k^-L f in[l + e(A'-fc")/^]-l[tan-i(:^--^) + tan-\'^-^)]dk' . (18) 
2n J ak' c c 

Integrating by part under the consideration of approximation that tan^^(A;/c) k/c for c 3> Vq, we have 

e=-A^' + fc^ (19) 

where 



Because of —^VQ—k?<^ 1 and c ^ Vo, we obtain up to the first order that 



2n{p + ph) = 1 + 7--!=^ + — • (21) 

L ^V^ - fc2 C 



From eq. (^,([l^) and (^l|), we obtain an analytic expression of pik): 



2.p{k) = [1 + \ J + ^][1 + e^-^'+^'y^r\ (22) 
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Obviously, the p(k) is a Fermi-like distribution. When c, L — > cxd, 

= i + ,(-!.+.^)/T - (23) 

which is just the distribution of free Fermi gas. The chemical potential /i should be positive-definite for the positive- 
definite density of roots. 

As all particles are bounded in the potential well, i.e. Max{k) ^ Vq, p{k) should vanish almost for fc > Vq . This 
requirement together with eq. (|2^ ) gives 

iV^ - p')/T » 1. 

The system being thermal equilibrium exhibits a simple dependence on the large momentum cut-off Tq (Tq — Vq — /i') , 
due to boundary effects. Therefore, the system can be in a state of thermal equilibrium only when T <C Tq. Otherwise 
some particles may overcome the potential energy at the boundary and escape out of the well. Substituting the 
obtained e into eq. ([T^), we obtain 

F^f,N~— [\n[l + e^^^'-''"'>/^] — [k+-sm-\ — )]dk. (24) 
27r J dk L Vq 

As fc <C Vo, we can replace sin"^(fc/Vb) by k/Vo, then 

2L 2 f^" k^ 



In the low-temperature condition, the free energy becomes 



where 

The /i' is regarded as a mandation of chemical potential according to eq. (p^). However, we are not able to get 
an explicit result for the specific heat by partial derivative of eq. because the chemical potential /x might be 

temperature dependent. In order to observe some properties of specific heat at low temperature, we let c ^ oo and 
let Ho denote chemical potential at zero temperature. The /io is determined by 



D = — 

27 



1 Z"^^" 1 2 

-/ (1 + - ^-^^ )dk. (28) 
Considering limr^o f-iT)/ p,Q = 1, we have 



2^2 



^ = MoA^--(l + y|r)(^A^f + ^)- (29) 
TT LVq 3 24fil/^ 

We find that the specific heat at low temperature is Fermi-liquid like 

ttT 2 

^v' = ^r^(l + 7^)T. (30) 
6fio ' LVo 

Thus the interaction between the particles plays an important role to their statistical properties though the system 
we considered is a boson system. This is a model belong to the class of Haldane's exclusion statistics p3|. 
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IV. THE WEAK COUPLING LIMIT 



Obviously, eq. (10) can be written as 

e= -fi + k^ - ^ I e-"l"le*'=" cos(fc'cj) ln(l + e-'''^)dk' dw. (31) 
27r J 

Because in weak coupling limit c <C 1 , we have 

e(fc) + e - Tln(l + e-"'^) - f{k, c), (32) 



where 



This leads to 



/(fc, c) = ^ y / Lil^fl^e^fc- cos(fc'^) ln(l + e(''"'=')/^)dfc'dc.. 

n=l 



Like-wisely, eq. (^) gives rise to 



e-./T ^ jg(-^+fc^-/(fe,c))/T _ (33) 



2^p;,(fc) = l + i— =|===+5(fc,c), (34) 



where 



P { A \Tt Th I I Th 

g{k, c)^J2 nl ^'^'^ cosik'Lo)dk'duj 



n=l 

With the help of ( p3| ) , we obtain the distribution function of roots 

Physically, it represents the distribution of quasi-niomenta of the system as a collection. Because the density of root 
should be positive-definite, ^ + f{k,c) must be always smaller than the corresponding fc^, particularly, ^ + f{0,c) < 0. 



So eq. (35) is boson-like distribution. If considering the boundary effects so that p{Vo) ~ we have Vq — ^ — /{Vq, c) ^ 
T. There is a large-momentum cut-off Tq ~ Vq — /(Vb, c) — such that the system can be in a thermal equilibrium 
only when T <C Tq . 

Now we consider the free energy, and only take account of the leading terms for small c, 

12 1 

27rp(fc) = [l + ]^ — -7,= . (36) 



The free energy (|l^) becomes 



F = ,N + T^J[l + \^J==^] - e^---'y-]dk. (37) 
Since fc ^ Vb, sin^^(fc/Vb) — k/Vo, we have 

As for free bosons (i.e., c — Q), the chemical potential should be nonpositive-definite and smaller than the energy 
of any particles. Moreover, for a fixed concentration {D — N/L), it is a function of temperature decreasing as the 
temperature increases according to eq. (p^). Hence when temperature goes to zero, the chemical potential will 
approach to zero from negative value. 
At low temperature, eq. (p8|) becomes 
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' .(e-/^), (39) 



where Lin{z) is polylogarithm function with Lz3/2(1) — C(3/2) and (^(x) is Riemann's Zeta function. Since Lin{z) 
can be expanded into series of z and the chemical potential /x should be zero at zero temperature, we neglect the 
Lz3/2(e''^"^)'s dependence on T. Then the entropy has the form 

S=^J^{l + ^)U,/,{e^n (40) 
4^7r LVq 

and the specific heat is proportional to T^/^ which exhibits bose-gas like behaviors. 



V. CONCLUSIONS AND REMARKS 



In the above, we discussed the thermodynamics of bosons in a one-dimensional potential well on the basis of the 
exact solution of the model. Using the strategy of Yang and Yang jl^, we studied the general thermodynamic 
properties of the system. We considered the problem in strong coupling limit and found that the behavior of the 
system at low temperature is Fermi liquid like even though it is a boson system. Therefore the interaction plays an 
important role. Meanwhile we obtained the specific heat which is linearly dependent on the temperature T. For the 
weak coupling limit, we found that the system behaves like free boson gas at low temperature. 
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